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Abstract. We study the weighted heat trace asymptotics of an operator of 
Laplace type with mixed boundary conditions where the weight function ex- 
hibits radial blowup. We give formulas for the first three boundary terms in 
the expansion in terms of geometrical data. 



1. Introduction 

An important issue for several decades has been to obtain explicitly the coef- 
ficients of the short-time asymptotic expansion of the heat kernel associated with 
a Laplace type operator on a m-dimensional Riemannian manifold M [15, 19]. In 
mathematics this interest stems in particular from the link between the spectrum 
of the operator and the underlying geometry of M [18], but it extends to basically 
all of Geometric Analysis [15]. In physics the heat kernel asymptotic expansion has 
been realized to be a particularly useful tool to determine various approximations 
of effective actions and the Casimir energy [2, 9, 13]. 

Instead of simply analyzing the integrated heat trace one often puts a weight in 
the evaluation of the trace, sometimes called the localizing or smearing function. 
This function is introduced for various reasons. First it allows one to obtain local 
information from the integrated one, therefore, most importantly it is possible to 
recover the local behavior near the boundary. Furthermore, it is this smeared 
coefficient that appears in the integration of conformal anomalies relevant for several 
physical applications, see, e.g., [6, 11, 19]. For smooth localizing functions the 
results for the first few heat kernel coefficients are available for several years now 
[16, 19]. A detailed analysis of what happens for singular weighing functions has only 
been started recently. In the context of the heat content asymptotics the weighing 
function plays the role of an initial temperature distribution. In the context of 
black hole physics singular conformal transformations play an important role when 
mapping black holes to their Penrose diagrams [5]. 

The heat content asymptotics of an operator of Laplace type with singular initial 
temperature distribution and with Dirichlet or Robin boundary conditions were 
investigated in [4]. A similar study of the heat trace asymptotics with singular 
weighting function and Dirichlet boundary conditions was performed in [3]. In 
this paper, we conclude this line of investigation by extending the results of [3] 
concerning heat trace asymptotics to Robin, and more generally, to mixed boundary 
conditions. We anticipate that also this singular setting will find its applications in 
physics. 

1.1. Operators of Laplace type. Let M be a compact Riemannian manifold of 
dimension m with smooth non-empty boundary dM . Let V he a, smooth vector 
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bundle over M and let D be an operator of Laplace type on the space of smooth 
sections C°°(y). This means that locally we may express D in the form 

(l.a) D = -{g^^d^^d,^ Id+A'^d,^ + B) 

for suitably chosen matrices A" and B where we adopt the Einstein convention and 
sum over repeated indices and where g^"^ denotes the inverse matrix. It is possible 
to express D invariantly [15] using a Bochner formalism. There exists a unique 
connection W onV and a unique endomorphism of so that 

where we use ';' to denote the components of multiple covariant differentiation. Let 
r be the Christoffel symbol of the Levi-Civita connection. We then have 

^ ' ' E = B- g^f^id^^cj, + Lj^cj^ - uj^T,^^) . 



1.2. Boundary conditions. We recall the formalism of Branson and Gilkey [7]. 
Let e > be the injectivity radius of the boundary dM in M. Use the geodesic 
flow defined by the unit inward normal vector field dr to define a diffeomorphism 
between the collar := dM x [0,e] and a neighborhood of the boundary in M 
which identifies dM x {0} with dM. The curves r — >■ {yo,r) for r S [0,e] are unit 
speed geodesies perpendicular to the boundary and r is the geodesic distance to the 
boundary. 

Let X G C°°{F,nd{V\aM)) satisfy = 1- Extend x to the collar so that 
Va^X = 0- Let Il± := ^(l ± x) be projections on the ±1 eigenbundles V± of x- 

Let S G End(y|aM) be an auxiliary endomorphism with II+S' = STI-i. = S. If 

(j) € C°°(V), let B — B{x, S) be the mixed boundary operator. 

(i.c) := {n_<^}|g^ © {n+(Va, + 5)n+^}|3^ . 

Let Di3 be the realization of D with this boundary condition. We set 11+ = to 
define the Dirichlet boundary operator Bxi and we set n_ = and 5 = to define 
the Neumann boundary operator Bj^f- 

Operators of this type arise when studying the Gauss-Bonnet theorem for man- 
ifolds with boundary [15] and will play an important role in the analysis of Section 
4. Let A = dS + 6d he the Laplace-Beltrami operator on the space of smooth 
differential forms. Let {y,r) be coordinates on the collar C^. Set 

n+(rfy/) = dyi and n_(dy/ A dr) = dyi A dr 

and define the absolute boundary operator Ba by taking 

Ba{4>id.y' + ipjdy' A dr} = {{dr<i>i)dy'}\aM ® {'4'jdy']\aM ■ 

Extend the second fundamental form L (see Section 1.6 below) to act as a derivation 
on the space of differential forms. Then: 

War{fidy') = {dr + L){fidy') so Ba=B{x,-L). 

Let Ap be Laplacian on the space of smooth p-forms. If M is a closed mani- 
fold, then ker(Ap^B^) is naturally isomorphic to the topological cohomology groups 
HP(M\ C). Relative boundary conditions Br are defined similarly using the Hodge 
* operator and one may identify ker(Ap_Ba) with the relative cohomology groups 
HP{M,dM;C). 
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1.3. The heat equation. For f > and (j) e Li^iV), let u = e~*'^^(j) be the solution 
of the heat equation: 

{dt + DB)u{x;t) = 0, Bu = 0, lim u{-,t) = <!){■) in L'^{V) . 

Let dvolM (resp. dvolaM) be the Riemannian measure on M (resp. dM). There 
is a smooth kernel PD,B{x,x;t) which gives the fundamental solution of the heat 
equation: 

u{x;t) = / pDB{x,x;t)(t){x)dvolM{x) . 

J M 

If is formally self-adjoint with respect to a fiber metric, we can take a spectral 
resolution {Aj^, f?,^} for Dg where {6^} is a complete orthonormal basis for L? with 
39^ = and D0^ = X^O^. We then have 

Pdb {x, x,t) = '^ e"*-^" 0^ {x)e^ (x) . 

1/ 

This series converges in the C°° topology for i > 0. (There are some additional 
notational complexities in the bundle valued case we suppress in the interests of 
simplicity). 

1.4. Weighting functions. We study the weighted heat trace Tri2 (Fe"*^''). Pre- 
vious work has concentrated on the smooth section - we review that work presently 
in Section 1.5. However, in this paper, we shall concentrate on a more general set- 
ting and consider the following class of smearing or weighting functions. Let a <1. 
Let F be a smooth function on the interior of M. We assume that r"F is smooth 
on the collar := dM x [0,e]; the parameter a controls the growth (if a > 0) or 
decay (if a < 0) of F near the boundary. We expand F in a modified Taylor series 
near the boundary: 

F(y,r) ~r-"(Fo(y) + rFi (y) + r2F2(?/) + ....) where 
Fi{y) = U9ry{r"F}\r=o. 

We remark that the assumption that a < 1 ensures that F e L^{M). With Dirichlet 

boundary conditions, the fundamental solution of the heat equation vanished to 
second order on the boundary and it was possible to consider the region a < 3; 
logarithmic singularities then appeared when a = 1,2. This is not possible in the 
more general situation since the fundamental solution of the heat equation pn^ need 
not vanish on dM and we must restrict to a < 1 to ensure convergence. 

1.5. Heat trace asymptotics in the smooth setting. Suppose a = so that 
F is smooth on all of M; this is the case considered classically. Work of Greiner 
[17] and of Seelcy [21] shows: 

Theorem 1.1. Let D be an operator of Laplace type on a compact Riemannian 
manifold M with smooth boundary. Let Ds be the realization of D with respect to 
the mixed boundary conditions B given in Equation (l.c). There is a full asymptotic 

series as t \.0 of the form: 

oo 

Trz,2(Fe-*^^) ~ (47r)-'"/2i-"/2 ^ i"a„(F, £)) 

oo 

+(47r)-'"/2t-(— ^ t'l^af{F, D, B) . 
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There are local invariants defined on M and on dM so that 

an{F,D)= / F{x)an{x,D)dvo\M{x), 
Jm 

I 

af{F, D,B) = Y, Fi{y)a\%y. D, B) dvoW(t/) • 

These invariants play an important role in index theory; they are also important 
in regularization results for mathematical physics [15, 19]. We remark that we have 

used a different indexing convention than is sometimes used in the literature and 
that we have handled the normalizing constants involving An slightly differently. 

1.6. Local formulas. One has explicit combinatorial formulas for these invariants 
in the smooth setting; the interior invariants are known for n < 5 [1, 14, 22] and the 
boundary invariants are known for ^ < 5 [8, 16, 19, 23]. We introduce the requisite 
notation as follows. 

Let Rijki be the components of the curvature tensor of the Riemannian manifold; 

with our sign convention, i?i22i = +1 on the unit sphere in R^. Let r := Rijji be 
the scalar curvature of the manifold. Let Pmm = Rimmi be the normal component 
of the Ricci tensor. Let Lab be the components of the second fundamental form 

on the boundary relative to an orthonormal frame {ei, em_i} for the tangent 
bundle of dM] Lab = gi^rm^ e^^b)- Relative to the coordinate frame, we have 

L{dfj,,d^) = rf,^rn = -^drg^^ for 1 < /x, < m - 1 . 

Thus ill = 1 for the unit disk in R^. Express D = - {g'^'Vv^^ + E) where V and 
E are as in Equation (l.b). Then: 

Theorem 1.2. 

(1) ao{x,F,D)=Tv{FlA}. 

(2) ai (x, F, D) = Ti{QFE + Ft Id} . 

Theorem 1.3. 

(1) a^''(2/,F,D,e) = ^Tr{Fox}, 

(2) a\^{y, F, D,B) = ^ Tr{2FoLaa Id +3F,x + 12^0^}, 

(3) al^iy, F, D,B) = ^ Tr{Fo[96xE + 16tx - SpmmX + LaaLbb{13U+ - 7n_) 

+LabLab{2n+ + 10n_) + 96LaaS + 192S^ - 12X;aX;a] 

+Fi[Laa(6n+ + 30n_) + 96^] + 48F2X}. 

1.7. The shifted asymptotic series. If a 7^ 0, there is a shift in the power of t for 
the boundary invariants but the interior scries discussed in Section 1.5 is unchanged. 
In [3] , we used the calculus of pscudo-difFcrcntial operators to establish the existence 
of an asymptotic scries with Dirichlct boundary conditions. The same approach 
extends directly to the situation at hand to yield the following generalization of 
Theorem 1.1: 

Theorem 1.4. Let D be an operator of Laplace type on a compact Riemannian 
m,a,nifold M with smooth boundary. Let Db be the realization of D with respect to 
the mixed boundary conditions B given in Equation (l.c). Let a < 1. Let F be 
smooth on the interior of M and let r"F be smooth near the boundary. There is a 
full asymptotic series as 1 10 of the form: 

00 

Tri2(Fe-*^«) ~ (47r)-™/2i-'"/2 ^ ra„(F, £>) 

n=0 

00 

+ (47r)-'"/2i-(--l)/2 ^ ^(l-a)/2^bd^^p^ jj^ _ 
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The interior invariants an{F,D) are as discussed in Theorem 1.1. There are local 
invariants a^'^ ^{y,D,B) which are real analytic in the parameter a, so 

r ^ 

JdM 

The interior invariants do not depend either upon a or upon B and are described 
for n = and n = 1 by Theorem 1.2 (and are known expUcitly in the hterature 
for n < 5); thus our attention will be concentrated on the boundary invariants and 
upon extending Theorem 1.3 to this more general setting. The analyticity of the 
invariants in a will play a crucial role. We shall often restrict to the case a ^ Z in 
proving certain identities to ensure that the interior terms and the boundary terms 
do not interact. We shall also often assume a << to avoid convergence questions. 
The result for general a will then follow by analytic continuation since the local 
invariants are real analytic in the parameter a. We set the local boundary heat 
trace density to be: 

e 

a\%y,F,D,B) := ^ F,(y)a,^;i,,(y, 

i=l 

1.8. Dirichlet Boundary Conditions. We computed the boundary invariants for 
Dirichlet boundary conditions in [3]; the following result is a consequence of those 
computations and forms an essential starting point for the study of the general case: 

Theorem 1.5. Let B define Dirichlet boundary conditions. Let '■— \T (-'-j^). 

(1) al%{y,F,D,B) = K„Tr{-JbId}. 

(2) aj^„(y, J^, A B) = Ka-iTr{-Fi Id + 2^^biaa Id}. 

(3) 4^„(y, F, D, B) = K,_2 Tr{-F2 Id+^^F,L,, Id 

+ qFoR 

amma 

~ 3(l-a)^oRijji Id - JZ^FqE} . 

1.9. Heat trace asymptotics for mixed boundary invariants. The follow- 
ing is the main result of this paper. It generalizes Theorem 1.5 to general mixed 
boundary conditions: 

Theorem 1.6. Let ag^a be the invariants o/ Theorem 1.4. 

(1) ag^„(y, F, D, B) = Tr{Fo[n+ - n_]}. 

(2) al'^Jy, F, D, B) = Tr{Fi[n+ - n_] 

+FoI'aa[ 2("_T)(3-a) ^+ + 2(^n_] + ^FqS} . 

(3) a|f„(y,F,Z),^) = fc„_2Tr{i^2[n+-n_] 

+Fi [Laa ( 2(a-~2^(4-a) ^4- + 2^^ II-) + 

+i^o[-gPmm(n+ -n_) 

,7- T ( a''-6a'^-a^-2a+104 tt (a-7) -pr ^ 

+^nn^'''A8(a-6)(a-4)(Q-2)(Q-l)^^+ ~ %{a-(>)^^-) 

,7- ,r r a-'^-10a^+21a+4 tt a-5 tt 1 

■^ljahi^ah{ 4(a_6)(a-4)(a-l)^^+ + 4(a-6) ^^"J 

2(a^-a 



(a-l)(Q-2)(a-4) ' (2-a)(l-a) 

+ 3(T^(^ + 6i^)(n+ - n-) - (,_i)\,_4) X;aX;a]}- 

Here is a brief outline to the paper. In Section 2, we express the invariants a^^, 
a^i"^, and a'2'a in terms of geometrical quantities with 8 undetermined universal 
coefficients ^J,; we refer to Lemma 2.8 for details. In Section 3, we determine the 
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coefficient of S in ai and the coefficient of in 02 by performing a computation on 

the interval. In Section 4, we examine absolute and relative boundary conditions 
in dimension 2 to derive additional relations. In Section 5, we use the calculus of 
pseudo-differential operators to complete the computation. 

2. The method of universal coefficients 

2.1. Weighted homogeneity and dimensional analysis. We assign weight k to 
the /e*^ derivative of the metric, weight A: + 1 to the A:*^ derivative of the connection 
form oj of Equation (l.b), and weight k+2 to the fc**^ derivative of the endomorphism 
E of Equation (l.b). We also assign weight k to the fc"^ tangential derivative of 
X, weight k to F^, and weight fc + 1 to the fc"^ tangential derivative of S. Thus, 
in particular, the components Rijki of the curvature tensor have weight 2 and the 
components Lab of the second fundamental form have weight 1. Standard argiiments 
using dimensional analysis shows establishes the following result; we omit details 
in the interests of brevity and instead refer to [3, 15, 19] where similar results were 
established: 

Lemma 2.1. The local invariants ^ of Theorem I.4 are weighted homogeneous 
of degree £ — i. 

2.2. Orthogonal invariants. Weyl's theory of orthogonal invariants [24] may be 
used to construct a spanning set for the space of invariants which are homogeneous 
of weight k. One uses the metric to contract indices in pairs. We let x-.a denote the 
components of tangential covariant differentiation of the tensor x- Lemma 2.2 then 
leads to the following result; again, we omit details as by now the arguments are 
standard: 

Lemma 2.2. There exist universal constants {^^^j so that: 

(1) al'^^iy,F,D,B) = TV{J^o[^°'+n+ + ^"'"n.]}. 

(2) a\<^^{y,F,D,B) = Tr{F,[gl;+n+ + ^i-H.] 

+Fo[Ka{ey[i+ + ^^'-n_) + g^^+S]}. 

(3) F, D, B) = Tr{F2[QyU+ + ^^.-n.] 

+F,[Laa{gl:+n+ + qI;-II_) + p^^+S] 

+i^o[Pmm(eI^+n+ + g^-U^) + LaaLbb{gl;+U+ + gl;-U^) 
+LabLab{g^^+Ti+ + e&-ii-) + E{gi°'+n+ + oi"'-n_) 
+T{gi^'+Il+ + gl'^-n^) + g]^'+LaaS + g^^S^ + gl^'+ X;aX;a]} ■ 

2.3. Product formulas. The following is a useful observation. 

Lemma 2.3. Let M = Mi x M2 where Mi is a closed Riemannian manifold and M2 
is a Riemannian manifold with smooth boundary.. Let Di be operators of Laplace 
type on Mi and let B be a mixed boundary operator on M2 which we extend to M. 
Then: 

ae'!a,i{i^i,y2),D,B)= ^ ak{xi,Di)a'^%^,{y2,D2,B) . 
2k+j=e 

Proof. Because the structures decouple, we have that e~*^^ = e~*^^e~*^^'^. Let 
F{xi,X2) = Fi{xi)F2{x2). Then: 

The desired result then follows by equating terms in the asymptotic series. □ 
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2.4. Dimension shifting. A-priori, the constants in Lemma 2.2 depend on the 
dimension. Fortunately, that is not the case. 

Lemma 2.4. The constants are independent of the dimension m. 

Proof. Let Mi and Di := —dg where 9 is the usual periodic parameter on 

the circle. Since Mi is a closed manifold, there are no boundary invariants. Since 
the structures are flat, a„(^, Di) = for n > 0. By Theorem 1.2, ao{6,Di) = 1. 
Consequently, by Lemma 2.3, 

It now follows that in dimension m is equal to in dimension m + 1. □ 

2.5. The coefficients of E and of r. In the proof of Lemma 2.4, we applied 
Lemma 2.3 with Mi = S^. We take a product with S"^ to establish: 

Lemma 2.5. We have the relations p^°'^ = p°'^ and p^^'^ = 

Proof. Let Mi = S'^ be the sphere of radius e in M'^, let A be the scalar Laplacian 
on S"^, and let Di := A — 5. We have t = and E = 5. We apply Theorem 1.2 
to see 

ao{x,D) = l and ai{x,D) = S + ^e~'^ . 

Let D2 be an operator of Laplace type on M2 and let B be mixed boundary condi- 
tions. We form M := Mi x M2 and D = Di + D2. We omit terms which involve 
neither 6 nor e and use Lemma 2.2 and Lemma 2.3 to see 

<a,0((a;i,2/2),AB) 

= Tr{S{gi''+U+ + gi''-U.) + 6-'{gi''+U+ + gi''-IL.) + ...} 

= Tr{(5+i£-2)(^o.+n+ + ^o.-n_) + ...}. 
Equating coefiicients of 6 and yields the desired identity. □ 

2.6. Degree shifting. In Lemma 2.5, we related p]^'^ and p^^'^ to p^'^. There 
are other relations of this form which are available: 

Lemma 2.6. 

(1) If'^<i< then ai^a,i = a^-i,a-i,i-i- 

jd = 



Proof. Choose x(^) to be a smooth cut-off function which is identically near r = e 
and which is identically 1 near r = 0. Let F^{y,r) := x{''')'''~°'i'^ f{y)\ we then have 
(-^a)j(y) = ^jf{y)- We suppress the interior terms to express: 



(2) 


91;^ 


o.± 

= 


(3) 




2.± 

= 


(4) 




0,± 

= 9a-2 



00 „ 

TVi2(i^^e-*^-) ^ (4^)-™/2i-(m-l)/2 y-^(£-a)/2 / j^^i^^ 



Let 1 < « < ^. Wc have = F^\. Equating powers of t in the asymptotic 
expansions then yields the relation of Assertion (1). We apply Assertion (1) with 
£ = 1 and i = 1 to derive Assertion (2); we apply Assertion (1) with 1 = 2 and 
z = 1 to derive Assertion (3); wc apply Assertion (1) with £ = 2 and i = 2 to see 
p^'^ = p^'^i and then apply Assertion (2) (after replacing a by a — 1) to establish 
Assertion (4). □ 
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2.7. Relating pure Neumann and pure Dirichlet boundary conditions. 

The following Lemma gives some relationships between the coefficients defining 
pure Neumann and pure Dirichlet boundary conditions. 

Lemma 2.7. ^0'+ + g^^- = g'^^+ + g];- = 0. 

Proof. Let M be the upper hemisphere of the unit sphere of M"*"'"^. Let 

be an isometric involution of 5"* whose fixed point set is the boundary of M. Let 

^(A, Asm) = {</.€ C^iS"") : Asn.(A = 

be the eigenspaces of the spherical Laplacian on 5™. Let {T*(f){x) := (j>{Tx). Since 
T is an isometry of S"", it commutes with the Laplacian and we may decompose 

E{X, Asm ) = E+ (A, As- ) ® E- (A, As- ) 

into the ±1 eigenvalues of T*. It is then immediate that the elements of E~ satisfy 
Dirichlet boundary conditions while the elements of E~^ satisfy Neumann boundary 
conditions. If pt> and p^T are the fundamental solutions of the heat equation of 
the Laplacians Ap and A^^^ on M for Dirichlet and Neumann boundary conditions, 
respectively, we may then conclude, after allowing for the renormalization of the 
norms of the eigenvectors, that: 

PD{x,x;t) + pj\j-{x,x;t) = 2ps'n{x,x;t) for xGM. 

Since 5™ is a homogeneous space, there are constants a„ so that 

ps™ {x,x;t)= psrr. (t) ~ (47rf) "^/^ ^ i"a„ . 

n 

Let F e C°°(5"") satisfy T*F = F. Then 

p OO 

Tri2{f e"*^^} + Tri2{f e"*^-^} ~ {4wt)-"'/^ / F{x) dyo\{x) ■ ^ Ta^ . 

We may suppose a ^ Z. Since there are no terms in the asymp- 

totic expansion of the right hand side of the above display, the boundary terms 
must vanish. The two relations of the Lemma now follow. □ 

The coefficients p^" may be evaluated using Theorem 1.5 After changing notation 
appropriately to simplify the relevant formulas, we summarize the results of this 
section in the following result: 

Lemma 2.8. There exist universal constants {1?^,} so that: 

(1) al'^Jy, F, D, B) = TV{Fo[n+ - Il_]}dy. 

(2) a\<^Jy, F, D, B) = Tr{Fi[n+ - n_] 

+FoL,aK^+ + 2^n_] + ^lFoS}dy. 

(3) a|^„(y, F, D, B) = K-2 TV{F2[n+ - n_] 

+Fi[L„„«_in+ + ^n_) + 

+Fo[-lpmm{^+ - n_) + i„„Lfe(,(^tn+ - i^n_) 
+ 4^n_] + i^lLaaS + 

+ 3(T^(^ + 6i^)(n+ - n-) + ^?|X;aX;a]}rf2/. 



HEAT TRACE ASYMPTOTICS 



9 



3. Special case computations on the interval 

We note that sr{s) = T{s + 1). Since «;„ := |r(i^), we have the following 
identities which we note for future reference: 

(3.a) = , = . 

Ka-2 1 — a Ka-l a 

Lemma 3.1. We have that i?^ = o,nd that ■d'l — -r^ — r. 

Prooj. Let M = [0, 1]. Let = r~" near a; = and i^a = near x = \. Let 
^ 6 e R. We form 

A:=a^ + 6, A*:=-a^ + 6, D := A* A = AA* = -{dl-\?) . 

Let D-D and D-ji be the realizations, respectively, of D with respect to Dirichlet 
boundary conditions and Robin boundary conditions with 5(0) = 6, and 5(1) = —h. 
Thus we may identify S-r,^ = A(^\qm- We integrate by parts to derive the Green's 
formula: 

\a*Au,v) - {u, A* Av)^{x)dx = | - {Au,v) + (u,^v)| . 

This vanishes if u = = on dM or if Au = Av = on dM. Consequently both 

D-p and D-ji are self-adjoint. If D(f> = 0, then (j)" = b'^<j) so <p = a^e''^ + aie~''^ . Thus 
(f> satisfies Dirichlet boundary conditions means (j> — and thus ker{Dx>) = {0}. Let 

{9^,X^} = {sm{TTiyx),2TT^iy^ + b^} for iy=l,2,... 

be a spectral resolution for Dx>. We have similarly that kei:{D-pi) = e~^^ ■ R and 
that {A*9v/\/\,,\i,} is a spectral resolution of D-j^ on ker(£>7^)-'-. Let and 
be the fundamental solutions of the heat equation for D-p and f-^, respectively. We 
compute: 



1/ 1/ 

dtpnix, x;t) = -J2 e-*^"A*eu ■ A* 6^ . 

This then yields the identity: 

29t{p2,(x,x;t) -pK(a;,x;t)} = -2^e-*^-{£)^, • - A*Q, ■ A* 9,} 

= 2 ^ e-*^" { {B'lB, - h'M,) + (C^^L - 26C^. + &^^.^.)} 

= dx{dx - 2b)px,{t;x,x) . 

We suppose a << and a ^ Z to ensure convergence and to ensure that the interior 
and the boundary terms do not interact; the general case then follows by analytic 
continuation. We integrate by parts to see: 

2at[Tri2{F„e-*^^} - Tri2{F,e-*^^}] 

= / 2Fo,dt{pv{x,x;t) - pTz{x,x;t))dvol{x) 

J M 

= I Fadxipx -2h)pT>{x,x;t)&vol{x) 

= [ {F'^ + 2bF^)pv{x, x; t) dvol{x) 
= TVi.{(F;' + 26i^;)e-*^-}. 
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Notice that d^Fa = —aF^+i and d^F^ = a{a + l)-Fc«+2 near the boundary of M. 
Since the underlying operator is the same, the difference of the interior terms cancel 

and we have: 

29i[Tri2{F„e-*^^} - Tri2{F„e-*^^}] 

~ - a)t(^— ^)/2{a^fjF«, I), B^) - 4%{F^, D, B^)} 
e 

= Tri2{(F4' + 26J'^)e-*^^} 

~ a{a+l)J2t^^~''~^y^a'^'!a+2iFo.+2,D,Bv) 

3 

-2a6^i('=-«-i)/24^+i(J'„+i,AB2,). 
fe 

We equate coefficients in the asymptotic expansions to see 

.3 - a){a\'^^{Fo,,D,B-D) - a\%{F^,D,Bn)} 

^ - ' =a{a + l)a\%^^{Fo,+2,D,Bv) - 2ahaW^+^{Fo,+uD,Bv) • 

We remark that although the argument is superficially similar to that used in 

Branson- Gilkey Lemma 3.2, the outcome is radically different owing to the ne- 
cessity to include the parameter a; in particular, there is no interaction between 
the interior and the boundary terms. 

We have E = —6^ and S = b. We take i = 1 in Equation (3.b) to see: 

{1 - a){a1UFa,D,BT,) - a1UFa,D,Bn)} 



{a + l)o5f„+2(^a+2,^,SB) - 2abao,a+i{Pa+i,D,Bv) . 



This leads to the identity: 

(1 - a){-Ka-i'&'i}b = 2a{Ka+i}b. 
We apply Equation (3. a) to see: 

2 —2a Ka+i —2a —2 4 



(2-a)Kc.-2&'|^^-4^ =«"&^ + 



1 — aKa_i 1 — a a 1 — a 
Finally, we take £ = 2 to see: 

(2 - a) {F„ ,D,Bt,)- {F^ ,D,Bn)} 
= a{a + l)a^^„+2(^a+2,£',Bi,) - 2aba>i'^^^,iF^+u D,Bn) ■ 
This gives rise to the identity: 

6 

'3(1 -1^ + 2]) 

We use Equation (3. a) to solve this identity for to see: 

^" = T^- 2-al-a "(" + ^^^= (l-aK2-a) - ° 

4. Absolute and relative boundary conditions 

We establish the following result by generalizing the 1-dimensional construction 
of Lemma 3.1 to the 2-dimensional setting. 

Lemma 4.1. We have that 

- 2(a-l)(3-a) " 

(2) (1 - an^l - 8^1} - 2a{a + l){<+2 + ^^+2} = ^ + ^^^^g^- 
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Proof. We modify an argument from McKean and Singer [20]. Let M be a Riemann 

surface. Let denote the Laplacian with absolute boundary conditions Ba as 
discussed above. Let A**^ := A° © and let A°<* := A^ Let {0^, be a spectral 
resolution of A^^ on ker(A^^)-L. Then 

is a spectral resolution of A"*^ on ker(A°'^)-'-. One has: 

2dt{pAi- ^'^ ~ P^i" ^' ^) J' 
= 2dt e-'^" {{0.,0.)-{{d + 5)e^,),{d + 6)e,)/K} 

i/:A„>0 

= -2 5^ e-*^^{X,{e,,e,)-{{d + 5)6,),{d + 5)e,)} 

= -2 ^ e-'^''{{is.''je.,e.)-{{d + 5)e,,{d + 6)e,)}. 

v:\u>0 

Now comes a crucial point. The Laplacian A^ decomposes as the direct sum of two 
scalar operators on A''^. Let ^ be a function. Then 

2Ao0 = ~26.ii, 

2{d + 6)0 = 2de = 29,16^ + 26';2e^ 

2{(Ao0, 9) - {{d + 6)9, {d + 6)9)} = -29,^9 - 29,i9,i = -{9^).ii 
= Ao{9,9). 

Next, let = ^^e^ A be a 2-form. We compute: 

2A2e = -29.ue^ A 

2{d + S)& = 2Se = ~29.ie^ + 29.2e\ 

2{(A2G, 9) - {{d + 6)Q, {d + 6)Q)} = -26^9 - 29,i9,i = Ao{9^) 

= Ao(e,e). 

Consequently we have 

-2 J2 e-'^-^{{A''J9^,9^)-{{d + 6)e^,{d + S)9^)} = - ^ e-'^^ Ao{{9^,9^)} . 
We suppose a << and a^'L. We may then integrate by parts to see: 

I 

^ 4f„+2(Ai^a, A-.Sa)*^^-"-^)/^ . 

Equating terms in the asymptotic expansion then yields 

(£ - a - 1) {af^^iF^, A^, B„) - af^^iF^, A^\ 

= -aX^^2{AF^,A-\B^). 

We specialize to the case M is the disk of radius 1 in R^. Introduce the usual 
coordinates {R,9) so that x = Rcos9 and y = Rsm9; the distance to the boundary 
is then given by r = 1 — i?. We have 

F„ = (l-i?)-", 

A = -d%- R-'dR - R-^dl 

AFa = -a{a + l)Fa+2 - aR-'^F^+i. 
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We have that = (1 — r)~^ = 1 + r + .... Since only the first 3 terms in the 
Taylor series expansion of AFa play a role for € = 1, 2, we obtain the identity: 

= aa,';i+2((a + l)F„+2 + Fa+l+F„,A'=^^?„) for £ = 1,2. 

We first set £ = 1. After canceling the factors of —a from both sides of the 
equation we get the relation 

The operator Ba on functions is the Neumann boundary operator and the operator 
Ba on 2-forms is the Dirichlet boundary operator. Near the boundary, wc decompose 
a smooth 1-form 6 = Qidr + Q2d6. The operator Ba on Qidr is the Dirichlet 
boundary operator and the operator Ba on 62^^ is the Robin boundary operator 
with S = —L. Substituting this into Lemma 2.8 yields the relation: 

„ , , f 1 a - 2 ,1 

We continue our computation: 

2{a + l)dl+, 



2 _ Aa + {a + l){a-2) _ + 3a-2 
1 — a 1 — a ' 

+ 3a - 2 



,91 

2(a + l)(l-a)' 



2(Q-l)(3-a) 2(a-l)(3-a)' 
Next we take £ = 2. This yields the relation: 

(1 - aM'^^{F^,A^^,Ba) - A°^SJ} 

Applying Lemma 2.8 with S = -1, Loa = LabLab = LaaLbb = 1 and X-.aX-.a = 8 
yields the relation: 

= a«„ + ^ + ia + l)(i?^+2 - + €+2 + 4^} • 

Since ^ = i^, K+i = and = (^^^r^^, we have 

(1 - am - (23^13^ - 8<} 

2« { 1^ + 2f^ + (" + 1) (<+2 - + ^^+2 + if^) } • 

This leads to the relation: 

(1 - afi^l - Ml} - 2a{a + l){<+2 + ^^+2} 

_ 8(a-l)-a^-a+4+a(a-3) , a(a+l)(a-l) 
~ a-2 4(a-4) 

_ 4(a-l) I a{a+l)(a-l) 
~ a-2 4(a-4) 

The desired result now follows. □ 

5. The pseudo-differential calculus 

In this section, we will use the pseudo-differential calculus to complete the cal- 
culation. Only the invariant X-.aX-.a genuinely involves a vector valued context; it 
will be determined by Lemma 4.1 once the remaining coefficients are determined. 
Thus we will restrict our attention to the case in which 

D = Am = -g'^^'d^d, + , 
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is the scalar Laplacian. We shall work with Robin boundary conditions. 

We begin by reviewing some fairly standard material. Let a = (ai, am) be a 
multi-index. We set: 



lal = ai + ... + am, a\ = ai! x 



m"*> C - (air) ^ - ^ (afr) 



a;" = x"^ X ... X x: 

We apologize in advance for the slight notational confusion involved with using a 
to control the growth of F and also to using a as a multi-index. We use the metric 
to raise and lower indices; "," will denote partial differentiation. 

We refer to [10, 12, 15, 17, 21] for additional material about pseudo-differential 
operators. We wish to construct the resolvent (Am — A)~^ for large A. We first 
suppose M is a closed manifold. In the evaluation of the heat equation asymptotics 
homogeneity properties of symbols are relevant and it turns out that collecting 
terms according to homogeneity is useful; the complex parameter A has weight 2. 
Expand the symbol of Am — A in the form a2{x, ^, A) -|- ai{x,£,) + ao{x,£,) where: 

ai{x,^,X) = \/^6''^^, and ao{x,^,X) = 0. 
We formally expand the symbol of the resolvent in an asymptotic series: 

oo 

(5.a) <7((Am - X)-^){x, ^,X)^J2 ^, A) . 

;=o 

The Qk are then determined by the recursive relations: 

1 = a2{x,^, X)q-2{x,^,X), 
(5.b) ^ l^[dlaj{x,^,X)][D%q_2-i{x,^,X)] for fc > 1. 

k=2+l+\S\-j 

To complete the proof of Theorem 1.6, we must examine q-2, and q^^. We 
summarize the facts we shall need and omit details in the interests of brevity - 
the fact that Am is scalar plays an essential role. Let greek indices range from 1 
through m. We have: 

q-2 = a2"\ 

9-3 = -019-2 +£-3,39-2: 

9-4 = -ao9-2 + C-4,39-2 + C-4,49l2 + C-4,59^2> 

where 

C_3,3 = -V^(5^"a2)(5>2), 

c_4,3 = aj ~V^id-^a,){d:a2) - (9^-a2)(a>i) - i (9^^^ aa) (9,^^02), 

C_4,4 = -3aiC_3,3 + V^{d^a2){d:c.3,3) + a2)(9>2)(5>2), 
3,3- 



o 2 

C_4 5 = 3c_ 
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One has that: 

q-4{x,C,X) = 

If the manifold has a boundary the expansion (5. a) has to be augmented by a 
boundary correction. To formulate the conditions to be satisfied by the boundary 
correction we expand about r = 0. One may express the metric on the collar in 
the form 

dsli = 9<Te{y, r)dy'^ o dy^ + dr^ . 

The coordinate y locally parametrizes the boundary, and r is the geodesic distance 

to the boundary, so 3; = (y, r). A tilde above any quantity will indicate that it is to 
be evaluated at the boundary, that is at r = 0. Furthermore, we use ^ = (w, t). 
We find 



°° Qk 
fe=0 ■ |5|<2 



i-=0 



with the notation 



D^r = n 



Introducing 
aj{y,r,w,Dr,X) = 



E|ai=2«5(y,r) ni=i 



m—l ui 



for i = 0,l, 
A for j = 2 



we define the partial symbol 

1 " d'^ 

a'(AM - A) = ^ pr'^ ^ (y, r, u;, D^, A) 

As it turns out, the symbols 



k\ 

fe=0 j=0 



r=0 



a«(2/,r,u;,Z)„A) = ^ 



i=0 k=0 

l-k=j 



1 . d'' 

— r 

A;! dr^ 



ai{y,r,uj,Dr,X) 



r=0 



have suitable homogeneity properties and using these symbols we write 

2 

a'{AM-X)= J2 a^'Hy,r,co,Dr,X). 

j=-oo 

We write the symbol of the resolvent as 
(t((Am - A)"^)(y, r, w, r, A) 



(5.C) 



= Yq-2-j{y,r,u},T,X) - e ^^'^^/i_2-j(2/,r-,a;,T,A), 

j=0 j=0 
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where the second term is the boundary correction. The factor e appears be- 

cause the operator constructed from these terms is the Op'{h) in [21], and Op'{h) = 
Op(/ie-^^''). This shows 

CSO 

o-'(Am - X) o'^h_2-j{y,r,co,T, A) = 0. 

3=0 

Here o denotes the symbol product on R™~^. Analogously to Equation (5.b) this 
equation leads to the differential equations 

= a^'^\y,r,u,Dr,\)h_2{y,r,uj,T,X), 
= a^'^\y,r,u:,Dr,X)h-2-j{y,r,L0,T,X) 

+ E 4f [£>>W(y,r,a;,£>.,A)l [(x/^D^)5/i_2-i(y,r,w,r,A) 



3,fc,i<j 

j=l+2+\a\-k 

For the present considerations we need h-2-j for j = 0,1,2, and we have more 
explicitly (repeated letters a, 6, c, ... run over tangential coordinates {1, 2, m— 1}) 

= a'-^'^{y,r,uj,Dr,X)h_2{y,r,uj,T,X), 

= a'^' {y, r, uj, Dr, X)h^3{y, r, uj, t, A) + a*^^' {y, r, ui, Dr, X)h^2{y, r, lo, t, A) 

+ [Dta<-^\y,r,o,,Dr,X)\ [iV^D'')hh^2iy,r,uj,T, X)] , 
= a^^^ {y, r,u,Dr, X)h-4{y, r, w, r. A) + o/"' (y, r, uj, Dr, X)h-2 {y, r, uj, t. A) 



+ 



{y, r, UJ, Dr, A)] [{y/^D'')bh-2{v, r, uj, t. A) 
Dt^a(2) {y, r, uj, Dr, A)l \{^/^D'')i,h-2{y, r, uj, t. A) 



+a^^^ {y, r, UJ, Dr, X)h-3{y, r, u), r. A) 



+ 



£»^a(2) {y, r, uj, Dr, A)l \{V^D'')bh-z{y, r, uj, t. A) 



The relevant equations for a^^^ {y, r, uj, Dr, A), z = 0, 1, 2 are 

a''^\y,r,uj,Dr,X) 

= a2{y,r,uj,Dr,X)\r=o = 9°'^'^a'^b + Dl - A, 

a«(y,r,w,D„A) 

= r{dra2{y, r, uj, Dr, A))|,.=o + 01(2;, r, uj, Dr, X)\r=o 
= rgfuJaUJb + \/^6"wa + yf-l¥Dr 

a^^\y,r,uj,Dr,X) 

= ^r^{d^a2{y, r, uj, Dr, X))\r=o + r{drai{y, r,uj, Dr, X))\r=o 

+ao{y,r,uj,Dr,X)\r=o 

= ^r^g^^aOJb + rV^b'^r'^a + r^/^b'^rDr■ 

The differential equations have to be augmented by a growth condition 
(5.d) h-2-j{y,r,uj,T, X) ^ as r — )• 00, 

and an initial condition corresponding to the Robin boundary condition 

B<j) = {dr + S)<p 
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considered here. The first few boundary symbols satisfy 

drh-2{y,r,UJ,T, X)\r=0 = V-^Tq-2{y,r,u>,T,X)\r=o, 

drh-3{y,r,uj,T,X)\r=o = -Sh-2{y,r,w,T, X)\r=0 + Sq-2{y,r,u,T, X)\r=o 

+y/^Tq-3{y, r, uj, t, X)\r=Q + drq-2{y, r, r, A)|r=o, 

drh^4:{y,r,uj,T,X)\r=o = -Sh-3{y,r,uj,T,X)\r=a + Sq-3{y,r,uj,T,X)\r=o 

(5.c) +^/^Tq-4{y, r, cj, r, A)|r=o + drq-?,{y,r,uj,T, A)|,.=o- 

Once the symbols /1-2-j have been determined, their contribution to the asymp- 
totics of the trace of the heat kernel follows from multiple integration. As before, 
we suppose r^F £ C°°{Cs). The contribution reads 



El— a— m I 



i=0 



VL{y,F,^M)dy 



dM 



with 



(5.f) 



r]^{y,F,AM) = Yl I j ds 



dre 



where 7 is anticlockwise enclosing the poles of in the lower half-plane. The 

integral with respect to s is the contour integral transforming the resolvent to the 
heat kernel. Note that from (5.c) the contribution to the heat kernel is minus the 
above. 

As will become clear in the following, with A = ^ |wp 4- 
of the type 



T s, we need integrals 



AJ(r2+A2)^ 



-rA 



The T integration can be done using 



I 



So 



dre 



pkljn 
'-ab... 



(v^)'^(-l)'+'=7r / 1 d 



2A dK 



i-i 



[A' 



fc— 1„— rAT 



(V 



■f 

Jo 



fe— 1„— rAn 



Performing the A-differentiation, different f-dependent functions would occur. It 
is therefore desirable to first perform the f-integration before performing the A- 
derivatives explicitly. This is achieved by noting that {z = A has to be put after 
the A differentiation has been performed) 

(\/^)*=(-l)"+'=+i7r 



rpkljn 



(n-l)! 

1 d 
2AdA 



A 



(n-l)! 

J__d_ 
2AdA 



/ doj dse^'^^^ 

z=A 

T{l + l-a) [ doj r dse^^'^^^ 

«=A 



(A + ^) 
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We can proceed in general by introducing numerical multipliers Cnki according to 



1 d 



71-1 



2KdhJ (A + 0)'+i-« 
The s-integration is then performed using 



^nkl 



^l+2n—k—a ' 



27r 

The final w-integrations follow from 



C{y) ^ 

by observing that 

/ duj WaiUJa^.-UJa^ 

In particular 



due 



-g UlaUlb _ 



1 V d 



d 



1 J dy°-^ dy'^' 



■C{y) 



I 



[ duj e-l-l' 



m— 1 y — 

71" = yg, 



'99ab, 



^ V 9 {QabQcd + QacQbd + Qadgbc) ■ 



Introducing the numerical multipliers dkijn according to 

2(^/=T)fc(-l)"+fc+l7^2^(^ + 1 - a)c„ki 



(n-l)!r( 



j-\-l — k—a 



+ n 



) 



we obtain the compact-looking answers 

jn / duj WaWfo-.-e"!"! 



rpkljn J 

^ab... ~ "'kljn 



where the last cj-intcgration is performed with the above results. 

Note that the numerical multipliers dkijn are easily determined using an algebraic 
computer program. Therefore, all appearing integrals can be very easily obtained. 

Let us apply this formalism explicitly to the leading orders, and we start with 
h-2{y,r,uj,T, X). The relevant difFcrcntial equation reads 

(a2-A2)/i_2(y,r,w,r,A) = 0, 

which has the general solution 

h-2{y, r, uj, T, A) ^ gie-"^ + g^e''^. 

The asymptotic condition (5.d) on the symbol as r — >■ oo imposes 92 = 0. The initial 
condition drh-2\r=o = \/— TT(j'-2|r=o gives g\ = — \/^(A(r^ + A^))^^r. Putting 
the information together we have obtained 



h-2{y,r,uj,T, A) 



-It 



-rA 



A(r2+A2) 

Performing the relevant integrals, with the notation 

[ dl= [ duj [ ds [ dfe^' (- f dre'^'' 
J Jk^-^ J-oo Jo \ Jj 
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produces 

dIh-2{y,r,uj,T, -V^s) = -V^dioiiTT' 



Taking into account the prefactor in (5.f) and the change of sign, this agrees with 
Assertion (1) of Theorem 1.6. 
In the next order we obtain 

- A2)/i_3(y, r, u>, T, X) = {E + [/i)e-'-^ + (F + U2)re-^'', 



where 



~ r2+A2' ~ A(r2 + A2) ' 

~ A(r2 + A2) + A3(r2 + A2) + A(V2 + A2)2 ' 



A2(t2+A2)' 



Note, for later arguments, that Ui{lu) and U2{lo) arc odd functions in co. Fur- 
thermore, for the scalar Laplacian at hand 6" = g'"^ri)c"'; thus they contain only 
tangential derivatives of the metric. 

Using for example the annihilator method, we write down the general form of 
the solution to this differential equation as 

h^s{y, r, w, r, A) cie^"^^ + C2re~''^ + csr^e^''''^ + c^e'''^. 

From the asymptotic condition (5.d) we conclude C4 = 0. From the initial condition 
given in Equation (5.e) we obtain 

= -^(^ + t^2)-^(F + C/l) 



A2(t2+A2) A(t2+A2) A(r2+A2)2 

6^2 2gfuj'=iOaUJbT 2gfu}aU}bT^ gfcJaCJb 



A(r2 + A2)2 A(t2+A2)3 A(r2+A2)3 A(r2+A2)2- 
Prom the differential equation we derive 

C2 = -^iF + U2)-^{E + U,), 

Collecting the available information, we see 

h-3{y, r, uj, T, A) = De"''^ + Sre"''^ + Cr^e"''^ + 0{uj), 

with 
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and where 0{uj) is an odd function in co. Furthermore, 0{uj) contains only tangen- 
tial derivatives of the metric. We next perform the multiple integrals; note, odd 
functions in w do not contribute. We obtain 

J dIh^3{y,f,uj,T,-V^s) 

= T^^^ V^g'^r9ab |--^^C^1021 - \d20il — ^^^^1111 + ^^^^1041 
— £^2013 + 5£^0012 + -^^^£^1131 + -^^^^1221 1 

-|-6'Vff7r^(-\/^rfio2i - rfooii) 

= Vg^^r (1 - t) ( 4(:!Tr(-l3) g°'ga6 + i^iS) . 

This confirms the value of in Lemma 4.1 and of i?^ in Lemma 3.1 after taking 
into account the prcfactor in (5.f) and the fact that g°'^gab = —g°'^gab,r = ^g'^^Lab- 
Up to this point the calculation can be considered a warm up for the next order. 
Leaving aside the 5- terms for the moment, we would like to determine the universal 
coefficients of the geometric invariants LaaLbb and LabLab- In terms of the metric 
these are determined by 

r _ _i~ 

^ab — 2^^^'^ ■ 

Using the Christoffel symbols 

^jk = Is*' i9lj,k + 9kl,j - 9jk,l) , 

and taking into account that with our sign convention the scalar curvature is given 
by the contraction g-'^Rijk' , we may expand the Riemann curvature tensor in the 
form: 

p I "p I p I j.'p ^"p ^ p '"p ^ 

^ijk — ^ jk ,i ik ,j ~t~ ^ in ^ jk jn ik 

The normal projection of the Riemann curvature tensor reads 

^amma — 29,ryo,c^r 2^ yac,rr 4^9 9 ,r ycaUbd 

2. n"-^ n'^'^ n n 2-n'^'^n 

— 4.y g yac.rybd^r 29 gac,rr- 

The above results suggest a strategy for the calculation. It suffices to consider the 
special case where the metric is independent of y. As a consequence, our answer 
will have the form 

(4^)-"/2 {Hr^-gac.rr + K -g'" 9ac.r~gbd.r + Lr'g'^gab.r'gcd.r} 

plus terms involving S. This has to be compared with the terms in a2^^{F,AM) 
that possibly contribute to these geometric invariants. In detail one can show these 
terms are (mod terms with tangential derivatives of the metric) 

3(1-0;) '^ ~ gPmm + l^t^aaLbb + l^aLabLab 

= ^12 ~ 3(l-c«)) 9'^'^9ac,rr + (~ 12(l-a) i^a) 9°'^9'^'^9ab,r9cd,r 
+ (4(1^ -M + ¥1) 9'''r''9ac,r9bd,r- 

So once we know H, K, L, as a check we can verify that 

^ = {t2~ 3(l-a)) '^«-2' 

and we can deduce 

'^a Ka-2 ^ 3(l-a) ' ^oc 1-a ^ 6- 

In summary, when writing down the differential equation for h-i{y,r,uj,T, X), we 
can neglect all terms that are odd in w as well as all terms that contain tangential 
derivatives of the metric. We obtain (up to irrelevant terms) 

{d^ - h?)h-^{y, r, oj, r, A) = Me"''^ + Nre'''^ + Pr^e'''^ + Qr^e"''^, 
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where 



M 
N 



-K¥D + ¥B, 



A2 



gfua^bD - K¥B + 2¥C, 



P = - 



2A(t2 + A2) 



Q 



gfuJa^^bC, 



with B,C,D given above. 

So the solution has the form, taking into account the asymptotic behavior (5.d), 

To simphfy the notation, let S := + A^. From the initial condition we obtain, up 
to irrelevant terms. 



a 



+- 



P_ SD_ 

A IT ^ AS2 ^ aP 

AS3 

6V^¥gfuJaUJbT^ - AV^gf^cOaiObT^ - 2V^gfg';fLOaUJbOOci^dT 



A54 



12y/^g-^'g-dujgOJbiOcCJdT'' 



+ 



AS5 ' AS2 

2y/^g°'^^LJaL0bT 6\^g'fg';fcJaUJbU)cUJdT 



AS3 



AS3 AS4 
From the differential equation we obtain the conditions 

M = -2A/3 + 27, N = -4A7 + 65, 
P = -6A6 + 12e, Q = -8eA . 

This determines the numerical multipliers /?, 7, 5 and e to be 



o__3_Q__l^_lJV_lM 
P ~ SA* 4 A3 4 A2 2 A ' 

A - -1 <3 - IP 

" ~ 4 A^ 6 A' 



^ - _3 Q. _ IP. _ IN 
I ~ 8 A3 4 A2 4 A ' 

- _i Q 



8 A- 



For Am, we have: 

¥'b 



1_ ~ab~ 
2.V yab,rj 

1 ;rab;rcd~ ^ 
4.9 .9 9ab,rgcd,r, 



gabg^r — 2-9 -9 9ah,rgcd,r, 

~ 2-9 -9 (Jcd,rgab,r " 2-9 f)ab,rr, 
9ab9,rr — ^9 9 9ab,rgcd,r - 5 gab,rr, 

gfg'^r {dabgcd + gachd + gadgbc) = g"^g'"''gab,rgcd,r + W^9'"^9ac,rgbd,r- 

Performing the integrations we obtain, modulo normalizing constants of 7r^'"~^^/^\/^ 
one obtains: 



oil 



' 9 9ah^rr 



-Id: 



'3013 



noi3 



+2\/— 1(^3014 — -^^Y^C^1012 



-Id 



1013 



+9°'^g'"^9ac,r9bd,T 



" ^^64^ ^1071 + 3C?2043 



5^0042 + -^^pC^lOSl 



8^^1031 — \/— 1(^3013 + \/— 1(^1013 — ^\/—ldiOli 

-\/^dioi4 + 6\/^d3oi5 + ^^dioi2 + 2\/^dioi3 - 3\/^dioi4 
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+g°-^g'"^gab,r9cd.,r 



— \d2mz + |c?0022 — -^^55^^1031 — 1^2022 + -^^^C^SOIS " -^^"J^rflOlS 
^^^^^3014 — ^^^-^dwii + 3\/— Ic^SOlS — -^^-T^C?1013 — ^''^ ^1014 



4^*^3014 — -^^-5-^C?1014 + 1'^3015 — -^^-J^C?1013 " ^^4^' 



— •^2022 — 2c?2023 + (^0022 + V^— 1^1012 + 2-\/^di013] 

+5*^ [—^—1^1031 — (^0021] I 



'^1141 + -^^-^^1121 



+9°-^g'"^gab,rgcd, 



14:1 I g -^IIZJ. 

— -^^^^1161 + 2C^2133 — 1^0132 + ■^^p^<ill41 — "^^^^1121 
— -^l^C^liei + g£^2133 — ]^C^0132 + -^^£^1141 + ]^ci2132 
— \d2ii3 + 1^0112 — -^^^ 



— "1231 



ii = g^-^gab 
+g°'^g'"^gac,rgbd, 

+ ^£^2223 — |c'0222] 
+9°'^g'"^gab,rgcd,r 



-dim + ■jc'oi3i ^ ^^j^dii2i — ^(ioiii 



^1211 



g^"1251 H g— "1231 §—"1211 32-"1251 



y yab,rycd,r J28-"1251 n 

+ 1^2223 — ^C^0222 — "^^^^1211 



-^1251 + ^3^cil231 — ■^^J^<il251 + ]^C^2222 



■^^£^1231 + 3C^0221 , 

I ~ab~cd~ ~ 
+ 9 9 9ac,r9hd,; 



■^^^c;i32i 



^■^^g^C^1341 + ■^^^'^1321 



~^i5?'^i34i + ^^d 



+SLa 
Si = .9"',9ab,r 

+g"'^g'^'^gab,rgcd,r 

£/ = g"^g'"^gac,rgbd,r §4- "1431 ~r y y yab,rycd,r J28""'1431 • 

Adding up all terms and simplifying using the functional equation and the doubling 
formula for the F-function, the contribution to the heat kernel coefficient reads 



-32-01321 

+ g°'^g'"^gab,Tgcd,, 



-^rfi431 



+ {T2-w=^)r'~9 

+ 49 g 9ab,r9cd 

This allows us conclude: 



i-o&acd;; ~ ( 1 _ 1 _ Q-'^-10a^+21a+4 

4,9 y 9ac,r9bd,r 6 4(Q-6)(a-4)(Q-l 

j_l^ab~cd~ ~ ( 1 I a^-6a^-a^-2a+104 A 1 
-|-4</ 9 9ab,r9cd,r y 3{l-a) ^ S{a-6){a-4.){a-2){a-l) J j ' 



"a 



a*-6a^-a^-2a+WA 
l(a-6)(a-4)(Q-2)(a-l) ' 

95 _ _ a''-10a'^+2lQ;+4 
a 4(a-6)(a-4)(a-l) ' 

,96 _ 2(a^-a-8) 

(a-l)(Q-2)(a-4) ■ 



(a-l)(Q-2)(a-4) ■ 

The values for a = reproduce the result for the smooth setting. We also confirm 
the result for i?^. We can now employ Lemma 4.1 (2) to determine d^. We find 

i98 — 1 

'^a — (a-l)(a-4) ' 

which in the limit a ^ reproduces the correct answer. 
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